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One of the greatest theoretical challenges in the build-up to the era of second-generation gravitational-wave 
detectors is the modeling of generic binary waveforms. We introduce an approximation that has the potential 
to significantly simplify this problem. We show that generic precessing-binary inspiral waveforms (covering a 
seven-dimensional space of intrinsic parameters) can be mapped to a two-dimensional space of non-precessing 
binaries, characterized by the mass ratio and a single effective total spin. The mapping consists of a time- 
dependent rotation of the waveforms into the quadrupole-aligned frame, and is extremely accurate (matches 
> 0.99 with parameter biases in the total spin of Ax < 0.04), even in the case of transitional precession. In 
addition, we demonstrate a simple method to construct hybrid post-Newtonian-numerical-relativity precessing- 
binary waveforms in the quadrupole-aligned frame, and provide evidence that our approximate mapping can be 
used all the way to the merger. Finally, based on these results, we outline a general proposal for the construction 
of generic waveform models, which will be the focus of future work. 
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I. INTRODUCTION 

The second-generation laser interferometric gravitational- 
wave detector Advanced LIGO is planned for first commis- 
sioning in 2014, and to reach design sensitivity in subsequent 
years [1-3]; Advanced Virgo J4j [5) and Kagra [6| are ex- 
pected to follow soon after. Current estimates of astrophysical 
event rates predict that the first direct detection of gravitational 
waves will occur in that time frame [7|. The coalescence of 
two black holes is among the strongest known gravitational- 
wave sources and a likely candidate for one of the first de- 
tections. The detection and subsequent analysis of gravita- 
tional waves relies strongly on the accuracy and completeness 
of theoretical waveform models. For black-hole binaries, this 
includes the inspiral, merger and ringdown of the final black 
hole, and current models combine information from analytic 
approximation methods and numerical-relativity (NR) simu- 
lations 0. 

To date, a number of theoretical inspiral-merger-ringdown 
(IMR) waveform models exist for non-spinning binaries and 
configurations where the spin angular momentum is either 
aligned or anti-aligned with the orbital angular momentum (a 
summary of these models is given in Ref. |8|). But most as- 
trophysical binary systems are expected to have arbitrary spin 
configurations, which lead to complicated precession effects. 
Although there does exist one preliminary precessing-binary 
IMR model the modeling of generic binaries remains a 
serious challenge. 

The complicated structure of precessing-binary waveforms 
suggests that in order to construct accurate IMR waveform 
models, we may need to produce numerical simulations that 
densely sample a seven-dimensional parameter space. At first 
glance, this does not seem feasible on the timescale of second- 
generation GW detectors (i.e., within the next ten years), al- 
though valiant efforts are underway 1 10 1. 

In this paper we introduce an approximation that has the 
potential to dramatically simplify the modeling of precessing- 



binary waveforms. Motivated by the results of our previ- 
ous work ifTTl . we show that the seven-dimensional space 
of intrinsic physical parameters of generic precessing-binary 
waveforms can be mapped to a two-dimensional space of 
non-precessing waveforms, parametrized by the mass ratio 
and an effective total spin parameter. The mapping consists 
of transforming the precessing-binary waveforms into a "co- 
processing" frame of reference, described by three rotation 
angles {y(f),/3(f),e(f)}. This is the "quadrupole-aligned" 
(QA) frame that we introduced in [11 1. The waveform mod- 
eling problem is then factorized into two much smaller prob- 
lems: (1) the construction of a non-precessing -binary model 
(and candidates for such a model already exist lfT2UT4 n. and 
(2) the construction of a model for the rotation angle func- 
tions {y(f), j3 (t ),e(f)} with respect to the binary's seven phys- 
ical parameters, which we expect can itself be further sim- 
plified. In this paper we do not address the (still very large) 
task of producing a model for the rotation angle functions, 
and we leave the behavior of the signal during the ringdown 
for future work. Here we restrict ourselves to an outline of 
the approximate mapping between precessing-binary and non- 
precessing-binary waveforms and test its validity on a series 
of inspiral waveforms generated by post-Newtonian (PN) the- 
ory. 

Generic binary systems undergoing quasi-circular inspiral 
are characterized by nine intrinsic physical parameters: the bi- 
nary's total mass M = m\ + /«2, the mass ratio q^m^frnx (we 
adopt the convention that m\ < ni2), and the six spin compo- 
nents Si and ^2- The total mass of the system sets the overall 
scale and can be factored out. The individual masses m\ and 
;«2 are uniquely determined given M and q. 

In the most simple cases, the two black holes are either 
not spinning, or their spin angular momenta are (anti)-aligned 
with the orbital angular momentum L. In these cases the in- 
spiral motion is confined to a fixed, time-independent plane. 
The orbital frequency of the motion then grows monotoni- 
cally with time and, when the GW signal is decomposed into 
spin-weighted spherical harmonics, most of the emitted grav- 
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itational energy is contained in the dominant harmonics, the 
(/ = 2, \m\ = 2)-modes. In the most general configurations, 
however, the two black-hole spins are not aligned with the or- 
bital angular momentum vector. Now the inspiral motion is no 
longer confined to a fixed plane. In these cases precession oc- 
curs, specifically two types: precession of the instantaneous 
orbital plane as well as precession of the spin vectors. The 
now far more complex motion is reflected in the emitted ra- 
diation in the form of strong amplitude modulations, which 
depend on the relative orientation of the binary towards the 
observer, and as a contribution to the binary's phase evolu- 
tion. These effects are illustrated further in Sec. Ill Al 

We have previously shown that precessing-binary wave- 
forms take a far simpler form when transformed into the 
quadrupole aligned (QA) frame ifTTI . In a nutshell, the QA 
frame approximately follows the instantaneous orbital plane 
of the binary. In this frame the binary is essentially viewed 
"face-on" throughout the course of its evolution. Note that 
this frame corresponds to a co-rotating, accelerated frame of 
reference. In this "co-precessing" frame the amplitudes of 
the waveform modes as well as their frequency evolution are 
significantly simplified and most of the energy is emitted in 
the (I — 2, \m\ = 2) modes, just as in a non-precessing bi- 
nary. In fact, in this accelerated frame the mode structure of 
a non-precessing binary appears to be restored (see Fig. 12 
in Ref. IfTTI for an NR example). It was this observation that 
suggested the idea that we pursue in this paper, that QA- and 
non-precessing-binary waveforms may agree well in both am- 
plitude and phase. Note that a related frame, defined by the 
direction of the Newtonian orbital angular momentum, was 
introduced in Ref. 1151 . along with the observation that the 
precession-induced phase oscillations can be removed in this 
"precessing frame". The key new result in this paper, beyond 
the use of the QA frame (which can be determined from the 
GW signal alone), is the simple mapping between precessing- 
frame waveforms and non-precessing-binary waveforms, as 
we describe below. 

In the context of gravitational-wave searches and param- 
eter estimation, waveforms from different binary configura- 
tions are most strongly characterized by their phase evolution, 
i.e., their rate of inspiral. When the black holes are widely 
separated their motion can be described well by PN methods, 
and in PN theory we see that the leading-order influence of 
the spin on the inspiral rate and the phase evolution is the 
spin-orbit coupling, S ■ L, where S is the total spin angular 
momentum of the two black holes, and L is the orbital an- 
gular momentum fl6l . If the binary precesses, then the pre- 
cession will introduce both secular and oscillatory changes in 
the phase, but in the QA frame, where the precession has been 
removed to some extent, we expect to recover the underly- 
ing orbital phase evolution, which will be similar to that for 
a non-precessing binary. Since the leading-order spin effects 
on the phase arise from the total black-hole spin, it is possible 
to make an approximate parametrization of non-precessing bi- 
naries by a single effective total spin parameter, XeS< an d this 
idea has been used in both inspiral [17| and IMR [ 1 3 1 mod- 
els. In this work we focus on inspiral PN models, and so we 
will use the same effective spin parameter as in Ref. IfTTI : see 



Sec. Ill For complete IMR waveforms other parameteriza- 



tions have been found to work better lfT2l[T3l . but in this work 
we restrict ourselves to PN inspiral waveforms. 

Our hypothesis, then, is that precessing-binary waveforms 
can be approximately mapped to non-precessing-binary wave- 
forms, and that the equivalent non-precessing-binary signal is 
parameterized by the mass ratio and Xeff- It is the goal of 
this paper to quantify the accuracy of that approximation. Our 
approach is to consider a selection of PN inspiral precessing- 
binary waveforms, and to match them against a family of non- 
precessing-binary signals, to determine the best-match value 
°f XeS- We can then see how well these values agree with our 
expectation, and the level of agreement with the best-match 
waveform. We use PN waveforms because they allow us to 
study the long inspiral regime, and they are far more compu- 
tationally convenient to produce than numerical simulations 
of only the last ~ 10 orbits before merger. 



In Sec. II A we will describe some of the general features of 
precessing-binary systems before giving a brief summary of 
the PN inspiral waveforms in Sec.|IIB| We provide a technical 



summary of the quadrupole-alignment procedure in Sec. II C 



We perform our study of PN waveforms in Sec. Ill and demon 



strate the efficacy of our approximate mapping not only for 
a large number of cases that exhibit simple precession, but 
also for an example configuration that undergoes transitional 
precession. In Sec. IV we show that the QA frame greatly 



simplifies the construction of hybrid PN-NR waveforms for 
precessing configurations, and discuss the potential extension 
of our approach to full IMR waveforms. Finally, we sketch 
a procedure to construct generic IMR models and discuss the 
issues that must first be overcome. 



II. PRECESSING-BINARY WAVEFORMS 

A. General features 

In this section we summarize the main features of 
precessing-binary systems and illustrate the effects of preces- 
sion on the gravitational-wave signal. For a comprehensive 
discussion of precessing-binary systems we refer the reader 
to Refs. lfT6l[T8l . which remain the standard references in the 
field. 

In non-spinning or spin-aligned cases the normal to the or- 
bital plane, i.e., the Newtonian orbital angular momentum Ljy, 
is well-defined and does not evolve and provides a useful di- 
rection with which to characterize the dynamics. In the pres- 
ence of precession any such characteristic direction becomes 
time-dependent. But one nearly-fixed direction in the binary 
configuration does exist: the direction of the total angular mo- 
mentum remains close to its limit when the binary has infi- 
nite separation. We denote this the "asymptotic total angular 
momentum direction" /_„ = J(t — > — °°). (We will use a hat 
to denote unit vectors). This is analogous to standard New- 
tonian solid-body mechanics, where the system rotates about 
the axis defined by the total angular momentum, which is a 
natural fixed direction. This is still true for Newtonian and 
first-order post-Newtonian binary systems. When spin effects 
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are included, starting at 1.5PN order, and in full General Rel- 
ativity, this natural direction of rotation still exists, but it is no 
longer fixed. The direction of the total angular momentum is 
now time-dependent. It evolves, but in cases with small pre- 
cession, and for large separations, it describes a precession 
cone that is rather small compared to all of the other time- 
dependent directions, like the orbital angular momentum or 
the spin directions fl9l . In a few special cases, where the 
orbital and spin angular momenta are nearly equal and oppo- 
site and the total angular momentum passes through zero dur- 
ing the inspiral, the direction of the total angular momentum 
changes rapidly; this is called transitional precession. 

The complex dynamics in precessing systems is reflected in 
the gravitational-wave signal. The precession introduces am- 
plitude modulations but also contributes to the phasing. Fur- 
thermore, the power contained in various gravitational-wave 
spherical harmonic modes, defined with respect to a fixed 
coordinate system, is significantly affected by precession, as 
power is transferred to the modes that were subdominant in 
the non-precessing configurations. Nonetheless, since the GW 
signal is to first approximation produced by the acceleration of 
the two bodies in orbit, the bulk of the energy is emitted along 
the direction of the orbital angular momentum L. This is the 
idea behind quadrupole-alignment: if we track the direction of 
the maximum energy emission, then we will also be tracking 
the orbital precession. 

This also provides insight into the signal observed from a 
fixed direction. The orbital angular momentum L precesses 
around 7_oo, and on average we expect the bulk power to be ra- 
diated in the direction of J^. This point is discussed in more 
detail in GUI . If a generic precessing binary now happens to 
be ideally oriented towards some static observer, i.e., the line- 
of-sight and coincide, only small amplitude modulations 
will be observed in the gravitational waveform, since the rela- 
tive orientation between the observer and the least-precessing 
axis of the binary does not change much. On the other hand, 
if the observer's orientation does not coincide with this most 
stable axis, then they will observe strong modulations. 

To illustrate this point, we consider a single-spin binary sys- 
tem with mass ratio 1:10, where the larger black hole has an 
initial spin of %2 = (0.75,0,0) and the smaller one is non- 
spinning. We expect to see only very few oscillations from 
directions close to /_<*>. This is illustrated in Fig. [T| the first 
panel shows the real part of the post-Newtonian (2,2)-mode 
of the GW strain h (see Sec. II B for details) of the precessing 
binary as seen along the direction of the binary's initial to- 
tal angular momentum. The second panel shows the real part 
of the (2,2)-mode for the same configuration, but for an ob- 
server whose line-of-sight does not coincide with the direction 
of the initial total angular momentum. In this case, the line- 
of-sight coincides with the direction of the Newtonian angular 
momentum at the beginning of the waveform. This direction 
varies with time and crosses the observer's line-of-sight af- 
ter each precession cycle. The amplitude peaks are observed 
when the maximal emission direction points towards the ob- 
server. (Note that since the line-of-sight coincides with L at 
t = 0, the amplitude has a maximum at this time.) 

The detectability of a signal in a matched-filter GW search 
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FIG. 1: The first panel shows the real part of the (£ = 2,m = 2) mode 
with J initially aligned with z, for the 1:10 binary described in the 
text. The second panel shows the same quantity, but now with L 
initial aligned with z. It is clear that if we are searching for signals 
with a monotonically increasing amplitude (as with non-precessing 
binaries), we may easily class the signal in the lower panel as a glitch 
in the data. 



can be estimated based on its best match against all of the 
theoretical signals in the template bank; the match is unity if 
the template bank contains precisely that signal, and is zero 
if the signal is orthogonal to all of the template-bank wave- 
forms. (Match calculations will be described in more detail 
in Sec. Ill ) For the two waveforms shown in Fig. [T] the first 
has a match against a non-precessing-binary template bank 
of over 0.97 (which means that more than 90% of such sig- 
nals would be found in a search), while the second has a best 
match below 0.7, meaning that it would most likely either be 
missed, or classed as a glitch in the data. These illustrations 
provide a useful perspective on the results in ifTTl . where it 
was shown that spin-aligned GW searches are very likely to 
miss a large number of generic signals at high mass ratios. We 
see that these signals are lost not simply because precessing- 
binary signals are very different from non-precessing-binary 
signals. Rather, for some orientations they look very similar, 
but the fraction of those orientations decreases as the preces- 
sion effects increase. 



These examples demonstrate the dramatic difference in the 
precessing-binary waveforms with respect to relative orienta- 
tion. But we have seen that even in the direction modula- 
tions remain in the waveform, and we will see in Sec.|IIC|that 
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these are further reduced when we go to the Q A-frame. 



B. Post-Newtonian Waveforms 

In order to produce the precessing post-Newtonian inspi- 
ral waveforms used in this analysis, we evolved the full PN 
equations of motion, which were integrated using a Mathe- 
matica package. In many studies of precessing binaries and 
in GW search work, it is common to use adiabatic inspiral 
models, for example the "SpinTaylor" equations, which are 
the precessing-binary extension of TaylorT4. But, in order 
to capture as much of the full physics as possible, we prefer 
to use instead an evolution of the full PN equations of mo- 
tion, formulated as the Hamiltonian equations of motion in 
the standard Taylor-expanded form, as we have done previ- 
ously I2TH231 . More specifically, we use the non-spinning 
3PN accurate Hamiltonian E4ti26l (see also G7H291 ) and the 
3.5PN accurate radiation flux [30-32]. We add both leading- 
order [[RO |33ti37l and next-to-leading order ||38H40l contri- 
butions to the spin-orbit and spin-spin Hamiltonians, and the 
spin-induced radiation flux terms as described in [41 1 (see 
also lfl6l [37l ). In addition we include the flux contribution 
due to the energy flowing into the black holes, which appears 
at the relative 2.5PN order, as derived in Ref. [42 1. 

The precessing PN waveforms were then generated mak- 
ing use of the explicit formulae for the waveform modes h( m 
as given in [43|. The positions, momenta and spins of the 
masses were read off the full PN solution and used to gener- 
ate the parameters for the construction of the precessing wave- 
form modes him- These modes contain only the leading-order 
spin contributions but higher-order corrections are contained 
in the dynamics, since the Hamiltonian is known to higher or- 
der (see above). Therefore, even if the h( m expressions were 
evaluated only at quadmpole-order, the waveforms would still 
show many features of precession, since the dominant contri- 
bution to the waveforms is from the motion itself. 

We have set up the source coordinate system as in Ref. |43 1, 
where Jq — (0,0, 1) and defines the total angular momentum 
direction at the initial separation. To achieve this, the PN ini- 
tial data {q,p, Si,Sz} were rotated by applying a standard ro- 
tation matrix about the y- and z-axes in the Cartesian source 
frame. This is purely a convention as all of the physics is in- 
variant with respect to rotations. The system was evolved for 
15M to reduce eccentricity (as done previously in numerical 
applications [21-23]), and then from an initial separation of 
Dj = 40M down to a final separation of Df = 6M. 

The orbital frequency of the motion is given by the general 
expression 
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where q is the relative separation of the point masses and q the 
relative velocity. The Newtonian orbital angular momentum is 
given by 



where [J. denotes the reduced mass 
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(2.3) 



The general PN orbital angular momentum vector L is given 
by 



L = q x p. 



(2.4) 



Note that Lpj and L differ significantly in the case of precession 
since q and p are no longer strictly parallel to each other, as 
explained in lfTTl[T6l unless the two masses m\ and mj_ are far 
apart. In the case of precession, the directions of L,Lm,Si and 
J are all time-dependent. The left panel of Fig. |6| shows the 
time evolution of J and L. We see that / evolves on a much 
smaller precession cone than the orbital angular momentum. 

The Newtonian orbital angular momentum Ln is defined by 
its polar coordinates {i(t),a(t)}, which are measured with 
respect to the z-axis of our Cartesian source system. The evo- 
lution of these two angles describes the dynamics of the in- 
stantaneous orbital plane. They are defined by 



(2.5) 
(2.6) 



The total phase of the binary is then constructed from the fol- 
lowing integral [43 1: 




4>(f) = / (a> orb (t')-a(t')cosi(t'))dt'. 



(2.7) 



The physical interpretation of the integral is as follows: the 
phase seen by an observer on the z-axis (which is the axis 
that defines our mode decomposition of the GW signal) is a 
combination of the actual motion of the orbital plane and its 
projection onto the xy-plane of the source frame. This can 
be better understood if we first consider a simplified example 
where 1 is constant (the orbital plane is tilted by a fixed angle 
with respect to z), and a is also constant (L^ precesses around 
/ with a constant frequency). Then we see that the average 
observed frequency of the objects' motion will be larger or 
smaller than (0 or b depending on the sign of a. It is the phase 
from this observed frequency that <t> describes. 

The symmetric and anti-symmetric spin combinations are 
constructed directly from the solution data: 



-Xh-fa) 



(2.8) 
(2.9) 



where the dimensionless spins %i ®& defined from the spin 
angular momenta Sj of each black hole by 



L N = n(qxq), 



(2.2) 
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Once all time-dependent dynamical parameters are con- 
structed, the waveform modes, h[ m , are evaluated. These are 
most commonly derived by expanding the complex polariza- 
tion h, 



h = k: 



(2.11) 



in the basis of spin-weighted spherical harmonics s Yi m with 
spin-weight s — —2 due to the nature of the gravitational field: 



Hm 



h{e,<p)- 2 Yt m {e,<p)da, 



(2.12) 



where * denotes complex conjugation. For a non-precessing 
binary this means that if the source frame was chosen such that 
L^v is parallel to z, the quadrupole contributions are I122 and 
hi -2- For precessing binaries, Ln is not in general parallel to 
z, and hence modes with m 7^ |2| appear even at quadrupole 
order. They only vanish when 1=0 and a = n. 
Schematically, the h\ m modes can be written as 



h\ m (t) = f(M,r,n,<a or b,l,<x,'P,Xs,Xa)- 



(2.13) 



The expressions are evaluated for a constant luminosity dis- 
tance r, which is scaled out of our results. Fig. [2] shows the 
magnitude of the (2,2)-mode for the same precessing case as 
described in Sec. [T] Despite this being a strongly precessing 
case (S-L = 0), long-timescale modulations are hardly notice- 
able. This is because a preferred frame was already chosen 
for the evolution, as described previously. Only an observer 
whose line-of-sight coincides with Jq will see a signal of this 
form. The appearance and strength of amplitude modulations 
strongly depends on the relative viewing angle, as illustrated 
in Sec. Ill Al 
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FIG. 2: The top panel shows the magnitude of the (2,2)-mode for 
the same strongly precessing case as in Fig.[T]over the whole length 
of the evolution, and over a length of the first 10000M in the lower 
panel. The source frame was chosen such that Jq ~ (0,0, 1), and L 
and X2 are initially orthogonal to each other with \\X2W = 0.75; the 
smaller black hole is not spinning. The close-up of the waveform 
magnitude over a shorter timescale reveals strong amplitude modu- 
lations. 



C. Quadrupole-Alignment 

The idea of quadrupole-alignment is to track the direction 
of the dominant radiation emission. This means that, to lead- 
ing order, it follows the precessing motion of the instanta- 
neous orbital plane. This allows us to significantly simplify 
the gravitational-wave signature by artificially removing the 
precession of the instantaneous orbital plane and describing 
the signal in a co-rotating way. In a previous study we have 
found evidence that the quadrupole-aligned direction actually 
tracks the full PN angular momentum direction, which differs 
slightly from the normal to the orbital plane. We will discuss 
this further in an upcoming publication l44l . 

We introduced the idea of the QA frame in ifTTI and il- 
lustrated its properties with reference to numerical-relativity 
waveforms. We also specified an explicit algorithm to de- 
termine the two time-dependent rotation angles {/3(f), 7(f)} 
that specify the direction that maximizes the amplitude of the 
[l = 2, \m\ = 2)-modes. A third angle, e(t), which adjusts the 
phase, was ignored in our original prescription, but its impor- 
tance was pointed out in |45|, particularly in whenever j3 is 
close to zero and y changes rapidly. An alternative algorithm 
to calculate these angles was later given in 11461 . 



The alignment itself is based on the general transformation 
behavior of spin-weighted spherical harmonics under coordi- 
nate transformations. This allows us to find the instantaneous, 
average direction of maximal emission by transforming the 
(£ = 2, \m\ = 2)-modes and averaging over the dominant har- 
monics. This direction is uniquely defined by two angles, /3 
and 7, which are determined by the maximization algorithm 
presented in IfTTI : 

(Pmax , Ymax) = max J\\kn(P,Y)\\ 2 + II hMPrfW 2 , 

< 2 - 14 > 

where h.22 and ^2,-2 are given by 

M/5,7) = t e- im '^d 2 m , 2 (-l5(t))h2„At), (2.15) 

m'=-2 

h-2(P,Y) = L e-^dlr^i-mhiAt), 

m'=—2 

(2.16) 

where d 2 , denote the Wigner d-matrices [47, 48 1. The 
maximization determines the two Euler angles Pmax and 
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Ymax- In general, the transformation of spin- weighted spheri- 
cal harmonics involves three degrees of freedom and, as noted 
in P31 , the third angle can be provided by the analog of 



Eq. (2.7 1, given the other two angles: 



e(t) 



Ymax (*')' cos p MAX (t')dt' . 



(2.17) 



We may set e(0) =0 without loss of generality. 

Once all three time-dependent angles (Pmax, Ymax, £) have 
been determined, the dominant quadrupole-aligned mode can 
then be written as 



-2ie(t) 



-im'rMAx(t 



4A- 



fiMAx{t))h 2m i{t). 



2 

I e 

m>=-2 

(2.18) 

All other QA modes can be constructed as well, as long as 
the /!/,„-modes for a given / are known. One may see that this 
transformation differs slightly from the one presented in ifTTI . 
This is because the numerical-relativity waveforms presented 
there are related to the PN waveforms in this work by an over- 
all complex conjugation. 

The three angles {y,/3,e} define a standard Euler rotation 
of the reference frame: a rotation by 7 about the z-axis, fol- 
lowed by a rotation by j3 about the >'-axis, followed by another 
rotation by e about the (new) z-axis. This is important to bear 
in mind if we consider the reverse procedure to "wrap up" a 
QA waveform back into its original precessing-binary form. 
In that case, the reverse procedure consists of applying the ro- 
tations in the opposite order, i.e., the same procedure but with 

{ 7 ,/3,£}^{-£,-j3,-7}. 

Although we expect QA waveforms to be useful tools in 
standardizing the representation of precessing waveforms for 
comparison purposes (as in, for example, Ref. [49 1 for equal- 
mass nonspinning waveforms), and in waveform modeling, 
they do not correspond to a signal seen by a gravitational- 
wave detector. The QA waveforms are the waves as seen in 
a very specific accelerated "co-precessing" frame. One of the 
consequences of this frame choice is that the usual relation- 
ship = —h no longer holds, as can be seen by inspection of 



Eq. (2. 18 1. Hence, in order to obtain quadrupole-aligned Weyl 



scalar modes, one has to construct the precessing modes first 
and then transform them into the quadrupole-aligned counter- 
parts. Note also that the QA angles will differ slightly when 
calculated from either h or ^4 (this point is also made in 11501 : 
the ¥4 angles tend to be smoother than the h angles). 

To leading PN order, the recovered angles correspond to 
the inverse Newtonian angles (i,a), but higher order contri- 
butions in the wave amplitudes lead to a deviation from those 
angles, which is consistent with the results from the pure nu- 
merical analysis in IfTTI . From that work we expect the angles 
we find to correspond to the smooth evolution of L in the limit 
of a complete description. The angles found by the maximiza- 
tion routine are shown in Fig. [3] They deviate slightly from 
the inverse Newtonian ones (—1, —a) due to higher-order PN 
contributions to the mode amplitudes. If we were to use only 
the quadrupole contribution of the h( m expressions, then we 
would recover the direction of Ljy. 

Once the three Euler angles are determined, those are 
then used to reconstruct the QA modes. Fig. [4] shows the 
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FIG. 3: The two panels show the angles found in an example of the 
maximisation routine. The first panel shows the inclination angle 
Pmax vs. time, the second panel shows the azimuth Ymax vs - time 
over the full length of the PN inspiral. 



quadrupole-aligned (2,2)-mode for the configuration shown 
in Figs. [T] and [2] 

In the next section we will present a detailed study of how 
these simplified QA waveforms compare with corresponding 
non-precessing cases. 



III. RESULTS 

The aim of this section is to test and quantify the accuracy 
of our hypothesis that generic inspiral signals can be mapped 
onto non-precessing counterparts (see Sec. [I}. Numerical- 
relativity waveforms are too short for a real inspiral com- 
parison and, moreover, it is computationally very expensive 
to produce a large number of accurate numerical precessing 
waveforms. Instead, we have restricted this analysis to PN 
waveforms to allow a more detailed study for a large subset of 
the precessing parameter space. 

First, we will take a look at simple precession and consider 
a range of spin configurations for two mass ratios. The first 
is q = 3 and includes the configuration of the numerical case 
that we studied in IfTTI . The second is q — 10, motivated by 
the observation that precession effects become more signif- 
icant for higher mass ratios; see, for example, Eqn. (2.11) 
in fl6l . and the results presented in lfl7l . We will show that 
the mapping works extremely well; the non-precessing wave- 
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FIG. 4: QA magnitude for the q = 10 configuration considered in 
Figs. [T] and [2] The top panel shows the complete waveform, while 
the lower panel zooms in on the first 10000M. We see that the oscil- 
lations in the amplitude have been reduced and simplified from those 
in Fig. [2] 



forms that agree best with each QA-transformed precessional 
configuration follow closely the XeS parameter that we dis- 
cussed in Sec. [I] (and will elucidate further below) and agree 
with them with matches above 0.99. Finally, as the most chal- 
lenging test of our hypothesis, we look at a case of transitional 
precession. 

This study covers only a small range of the full precessing- 
binary parameter space, but the configurations were carefully 
chosen to test the hypothesis for varying spin magnitudes and 
for two mass ratios within the range that is likely to be treated 
in IMR models in the near future, i.e., cases which can also be 
realized in current numerical simulations to high accuracy. 

From the PN expressions for the phase evolution of the bi- 
nary fl6l . we see that the dominant spin contribution is pro- 
portional to (S-L). We define K to be the angle between S and 
L: 



SL 



(3.1) 



When the spin interaction is restricted to the leading order 
spin-orbit coupling and radiation reaction is switched off, K 
is conserved and is a constant of the motion [ 16 1. When radi- 
ation reaction is included and, to a lesser degree, when higher 
order spin interactions are included, K has been observed to 
show only small variation in time. 



The agreement or disagreement between two waveforms 
is mainly due to their phasing. If the inspiral rate is signif- 
icantly different, two waveforms are not expected to agree 
very well. For the QA waveforms, the precession of the or- 
bital plane has been factored out, but the physical spins are of 
course present and contribute to the phase evolution. Thus, 
in general, we expect the best comparison waveform to be 
from a spinning-black-hole binary. At leading PN spin-order, 
where only the leading order spin-orbit term contributes, the 
spin-contribution is proportional to cos fc, and thus by look- 
ing at the leading-order terms, we expect that all waveforms 
with cos K = map onto non-spinning counterparts, while 
all waveforms with cos K ^ map onto spinning waveforms, 
which can be parameterized by an effective total-spin param- 
eter. This leading-order term can be represented by a single 
reduced spin parameter ifTTl : 



XeS — Xs. 



(mi - mi) 



%u,z 



76tj 

Tl3 



Xs 



(3.2) 



where 77 is the symmetric mass ratio. Note that this param- 
eter is not the same effective spin parameter as introduced in 
Ref. J36). In this work the effective total spin used is the re- 
duced spin parameter as defined by Eq. \i.2\ . 

In our study the non-precessing-binary comparison modes 
were parameterized by %i = Xi = (0,0,^). For each of these 
cases we have XeS = #(1 ~ 76tj/113). 

The first set of configurations was chosen such that K = 0, 
yielding an effective spin of zero. The second set was chosen 
such that all configurations have the same theoretical effec- 
tive spin of Xeff = 0.5, but with varying K. The details are 
listed in Tab. [I] and Tab. [TTJ The PN comparison family with 
(anti-)aligned spins was generated by the same method as the 
precessing ones, solving the full PN equations of motion and 
using the same h( m expressions |43|, where a = K and 1 = 0. 
This ensures that the results are not contaminated by differ- 
ences due to the choice of the PN approximant. 

The agreement between two waveforms can be quantified 
by a single number, the match which corresponds to a 
noise- weighted inner product (overlap) between them OP . 
Since QA waveforms are not in an inertial (detector) frame, 
and we are interested in quantifying the difference between 
two waveforms independently of a detector, we primarily use 
the white-noise spectrum S„(f) = 1. Match calculations are 
performed in the frequency domain and hence the FFTs of the 
time-domain waveform modes have to be computed first. The 
best match between two frequency-domain waveforms h\ (f) 
and 122(f) is defined as their normalized inner product maxi- 
mized over time and phase shifts (Af and A0): 



(hi\h 2 ) 

- max — 

aim y^TMMM 



where the inner product is defined by 



(hi\h 2 ) :=4Re 



'./in 



.An 



Hf)h* 2 (f) 

Sn(f) 



df 



(3.3) 



(3.4) 



In our examples the PN waveforms are defined in the fre- 
quency range fM £ [0.0018,0.01]. The upper frequency cor- 
responds to Mco sa 0.06, which is typical of the frequency at 
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which we would start using NR results in full IMR hybrids; 
in this study we are not interested in the performance of the 
PN waveforms beyond that frequency. Since the matches are 
calculated with a fiat noise spectrum, they are independent of 
the binary's mass. 

Although the QA waveforms are not in a detector's frame 
of reference, it is also instructive to calculate matches with 
respect to realistic detector noise curves. In this case differ- 
ent choices of binary mass correspond to giving extra weight 
to different frequency ranges in the waveforms, and provide 
a more stringent test on the robustness of our results. We 
repeated the match calculation for every configuration with 
the early Advanced LIGO ll52l and the zero-detuned high- 
power [53 1 noise curves. The matches were calculated for 
masses between 20M, S and 50M e in the frequency range be- 
tween 20 Hz and 8 kHz. 

The idea of the comparison is to find the non-precessing 
waveform as a function of % that gives the best match with 
each QA waveform of our study. If our hypothesis holds, then 
the best-match spin j;bm will be close to the effective spin Xeff- 
We present our results in the following subsections. 



A. Simple Precession 

The first two sets of PN configurations are cases of sim- 
ple precession. For most arbitrary binary configurations, sim- 
ple precession will occur and only a small set of configura- 
tions will undergo "transitional precession", as it requires fine- 
tuned physical parameters (see [18] and Sec. IIIB below). In 
the case of simple precession, the total spin angular momen- 
tum S precesses around the orbital angular momentum vector 
L and both of these vectors precess around the centre of the 
rather small precession cone described by Jq. This is illus- 
trated in the left panel of Fig. [6] 

Each precessing time-domain waveform was generated 
with respect to a source frame where Jo is initially parallel 
to the z-axis. The quadrupole-alignment algorithm was then 
applied to determine the time series of the two Euler rotation 
angles {f$MAx(t) ,YMAx(t)} ■ Given those, the third angle, e(f), 
was determined and Eq. ( |2.17[ i applied to reconstruct the time- 
domain quadrupole-aligned (2,2)-mode. 1 

The first set of configurations tests the mapping hypothesis 
for a vanishing proposed theoretical effective spin Xeff = 0, for 
various spin configurations for the two mass ratios 1 : 3 and 1 : 
10. The results in Tab.[I]suggest that the hypothesis works very 
well for single-spin systems with only the smaller black hole 
spinning. In these cases, we obtain best matches ^ 0.99 for 
the theoretical j e ff-value for both mass ratios. In the reversed 
cases, i.e., now the larger black hole is spinning, the maximal 
matches are still ^ 0.99 but we see a small parameter bias of 
Ax = 0.02. If both black holes are spinning with the same spin 
magnitude and the spins initially parallel to each other, the 



1 Higher modes can be reconstructed as well but here we consider only the 
dominant harmonic in the match calculations. 



parameter bias increases slightly to Ax = 0.03. Note that in 
all of these cases the match has a sharp peak at its maximum, 
but the match at the theoretical Xeff value is well above 0.97 
in many cases. 

The results do not change appreciably when the calcula- 
tions are repeated with the Advanced LIGO noise curves. The 
matches improve slightly as the mass is increased, but so does 
the bias in Xeff- However, the bias never increases by more 
than Ax = 0.01. The results for the 2OM bin are displayed in 
the last two columns of Tab. [T]and|ll] We would like to empha- 
size again that QA waveforms are not in a detector frame: the 
matches using the detector noise curves are only to rule out 
the possibility of spurious results with the white-noise curve. 

The second set was chosen such that all configurations have 
the same theoretical Xeff value, but that the amount of preces- 
sion changes due to a varying K angle. All configurations in 
this set are equal-spinning, i.e., the spins are initially equal in 
magnitude and parallel to each other. The results are given in 
Tab. [II] We see for both mass ratios q = 3 and q = 10 that the 
best-match x agrees with Xeff for small K. A bias appears as K 
increases beyond 30°, but is again never more than Ax = 0.02. 

The approximation that K is constant becomes less accu- 
rate as the binary approaches merger. Remarkably, the effec- 
tive spin value associated with the initial K value seems to 
characterize the best-match non-precessing-binary system in 
all cases. Even when using detector noise curves and choos- 
ing masses such that the late inspiral (when K changes fastest) 
is in the most sensitive part of the detector band, the best- 
match Xeff varies by only Ax < 0.04 from the value predicted 
by our hypothesis. However, it is likely that when we move 
to full IMR configurations, some other appropriate effective 
total spin will be more appropriate, as was found for the full 
IMR waveforms in Ref. [ 12 1. 

When interpreting these results, one should bear in mind 
that the phasing of a PN waveform can change significantly 
with respect to the choice of PN approximant. The matches 
that we calculated between QA and non-precessing wave- 
forms are in general far better than those between, for ex- 
ample, the same non-precessing configuration produced with 
TaylorTl and TaylorT4; see Fig. 6 in Ref. l54l . In this sense, 
our approximation can be considered to hold, well within the 
level of accuracy of our PN waveforms. 

We also emphasize once again that the QA waveforms do 
not correspond to the waveforms as seen by a detector, since 
the QA frame is accelerating, and would not be directly em- 
ployed in a GW search; the matches as shown therefore do not 
constitute a study of the efficacy of these waveforms for either 
searches or parameter estimation. What they do tell us, how- 
ever, is that if we were to take the non-precessing waveforms 
used in this study, and to apply the reverse QA procedure to 
them, i.e., "wrap them up" into mock precessing waveforms 
using the inverse QA angles calculated for each of these con- 
figurations, then we would expect them to agree well with the 
original precessing-binary waveforms. This study also sug- 
gests that if we were to construct a waveform model from 
"wrapped up" non-precessing waveforms, then it is possible 
that this model could be used to measure the effective total 
spin Xeff with only a small bias. However, the true behavior 
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q 


h 


Xi 


Xbm 




{Xbm) early 


(XBM)zdethp 


3 


(0,0,0) 


(0.75,0,0) 


0.02 


0.9815 


0.02 


0.02 


3 


(0.75,0,0) 


(0,0,0) 


0.00 


0.9997 


0.00 


0.00 


3 


(0.75,0,0) 


(0.75,0,0) 


0.03 


0.9576 


0.04 


0.03 


10 


(0,0,0) 


(0.75,0,0) 


0.03 


0.8209 


0.03 


0.03 


10 


(0.75,0,0) 


(0,0,0) 


0.00 


0.9999 


0.00 


0.00 


10 


(0.75,0,0) 


(0.75,0,0) 


0.03 


0.8075 


0.03 


0.03 



TABLE I: PN configurations with constant K = 90° and varying spins. The best matches, not necessarily for the predicted Xeff = but for 
the values displayed in column 4, are all well above 0.999 for q = 3 and above 0.995 for q = 10. ./#o denotes the match with the counterpart 
waveform that has Xefi — 0- The last two columns show the best match for two potential Advanced LIGO noise curves, evaluated for a 20M C ; ; 
binary. For all cases the best match is above 0.999 for both detector noise curves. 
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0.50 
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(0.444,0,0.529) 
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TABLE II: PN configurations with the same effective spin value Xeff — 0-5 but varying K for the two mass ratios 1 : 3 and 1 : 10. Xbm denotes 
the effective ^ e ff-value yielding the best match. In all cases the best matches are above 0.999 for q = 3 and above 0.997 for q = 10. ^#0.5 
denotes the match with the counterpart waveform that has Xeff = 0.5. Column 5 lists the match for the predicted ;£ e ff-value. The last two 
columns show the best match for two potential Advanced LIGO noise curves, evaluated for a 20Mq binary. 
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FIG. 5: The absolute value of the GW strain for a precessing binary, 
as viewed at an arbitrary inclination of 2.8 rad from Jo. The signal 
includes all I = 2 modes. The true signal (black) has the finer struc- 
ture; the other signal (red) was generated by twisting a non-spinning 
q = 3 waveform with the inverse QA angles. The dotted line shows 
the amplitude of the original non-spinning waveform. 



of such a model in a parameter estimation exercise requires an 
exhaustive study that is beyond the scope of this paper. 
To back up this claim, we performed the following exercise: 



from the first case in Tab. [I] we took the Xeff waveform, which 
is non-spinning q = 3, and wrapped it up with the reverse QA 
angles that we calculated for the {q = 3,Xi — ®>X2.z = 0.75} 
configuration. The resulting waveform is shown in Fig. [5j 
we have plotted the absolute value of the GW strain, con- 
structed from all I = 2 modes, at an arbitrarily chosen incli- 
nation of 2.8 rad from the initial direction of the total angu- 
lar momentum, Jq. Also shown is the same quantity for the 
"true" precessing-binary waveform, and for comparison we 
also show the original non-precessing-binary waveform, con- 
structed from only the (£ — 2, \m\ = 2)-modes. We see that 
the twisted-up non-precessing-binary waveform captures the 
main features of the amplitude of the true precessing-binary 
waveform extremely well; how well the phases agree can 
be judged by calculating the match between the two wave- 
forms. This we did, once again over the frequency range of 
fM <E [0.0018,0.01]. Note that now we are considering wave- 
forms as they would be observed in a detector. 

We find that the match between the true precessing-binary 
waveform and the mock-precession waveform have a match 
greater than 0.97 for all masses and binary orientations. By 
contrast, the match between the unmodified non-precessing 
q = 3 waveform and the true precessing waveform is below 
0.97 even for the best-performing orientation. These results 
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provide an important cross-check that we can indeed mimic 
the original PN precessing-binary signal by suitably trans- 
forming the signal from a non-precessing binary. 

As an aside, note that there is one mode of the precessing- 
binary signal that we cannot fully model in this way, the 
(I = 2,m = 0)-mode. In the non-precessing waveforms, the 
(2,2)- and (2,— 2)-modes are complex conjugates of each 
other. When this is true, the transformed (2,0) -mode will al- 
ways be real. This can be seen from inspection of Eq. ( |2.18 1. 
But in the true precessing-binary waveform the (2,0)-mode 
has real and imaginary parts; it is straightforward to produce 
an example to illustrate this from Sec. IV of Ref. B31 . In order 
to capture these effects, we would need to break the symmetry 
between the non-precessing h( m modes, which would require 
that our model include unequal spins — this is therefore one 
limitation of a single-effective-spin model. In practice, how- 
ever, the relative signal power in the (2,0)-mode will always 
be small, and we expect the other errors in this approximate 
waveform, for example in the phasing, will be more signifi- 
cant in practice. 



B. Transitional Precession 

In the previous section we have seen that our mapping 
works extremely well in cases of simple precession; in fact 
it can be considered to be an exact mapping within the error 
bars of the PN phasing. In this section, we demonstrate that 
it also works in the more extreme case of transitional preces- 
sion ||T8l . This second type of precession occurs when L and S 
are almost opposite and equal in magnitude and so \ J\ is small. 
During the inspiral, the magnitude of S hardly changes but 
since orbital angular momentum is radiated away, the magni- 
tude of L decreases with time. With the appropriate choice 
of parameters, the total angular momentum / is initially small 
and positive, but due to the loss of orbital angular momentum, 
decreases until it crosses the xy-plane of the Cartesian source 
frame, where it changes sign. See Ref. |18] for an extensive 
discussion of transitional precession. 

As opposed to simple precession, where Jq represents the 
least evolving axis in the binary's geometry, this direction 
changes significantly during the transitional phase, as shown 
in Fig. [6] In order to test the validity of our precessing— ^non- 
precessing mapping for a transitional-precession case, we 
have chosen one specific configuration with PN parameters 
q = 10, initial separation D, = 53M, and initial spins %\ = 
(0,0,0) and = 0.65 • (0,-sin(3°),-cos(3°)). This is a 
single-spin configuration, where the initial spin is 3° from 
complete anti-alignment and the generated inspiral waveform 
is about 2 • 10 6 M long, terminating at a final separation of 
Df = 6M. 

It is worth mentioning that in order to produce a transi- 
tional phase, the parameters have to be fine-tuned such that 
J changes sign. If S and L were completely anti-aligned, no 
precession would occur at all. The transitional phase is not 
brief: it takes up most of the duration of the inspiral that we 
have calculated, and, as noted in Ref. lfT8l . cases where a bi- 
nary undergoes transitional precession within the sensitivity 



band of ground-based detectors are expected to be very rare. 

The dramatic change of the direction of J is reflected in 
the GW signal and the transitional waveforms in the standard 
source frame look particularly distorted when the total angular 
momentum crosses the xy-plane, as is shown in Fig. [7] 

We do not expect any of these features to be present in the 
quadrupole-aligned waveform, since we now track the direc- 
tion of dominant emission, and this is completely independent 
from any asymptotic direction of /. We see in Fig. [7] that this 
is indeed the case. The angles found by the maximization rou- 
tine are shown in Fig. [8] The zero-crossing of the total angular 
momentum occurs at t = 1.587 ■ 10 6 M, which is indicated in 
the figures with a vertical line. 

If our hypothesis were correct, then the QA waveform 
would be very close to a non-precessing waveform with j e ff = 
—0.572, from Eq. (3.2 1. As before, we compared the QA 



mode with a series of non-precessing waveforms with varying 
spin parameter to locate the non-precessing configuration that 
agrees best with the QA waveform. We find the best match 
to be 0.998 for a spin anti-aligned waveform with effective 
spin parameter %eff = —0.576. This is remarkably close to the 
theoretically expected value, with a bias of only = 0.004! 

On the other hand, naively using the non-aligned 
transitional-precession waveform and calculating the matches 
with the same comparison waveforms gives the same effective 
spin value, since the phase is dominated by the inspiral rate, 
but yields a best match of only 0.940. Note also that this is for 
the (2,2)-mode as seen from only one orientation; for many 
other orientations that matches are likely to be far worse. 

This example demonstrates that even in the case of transi- 
tional precession, our method proves to be accurate (expected 
Jeff-value) and robust > 0.99) for mapping precessing 
waveforms onto single-spin-parameterized non-precessing- 
binary waveforms. 



IV. PN-NR HYBRID WAVEFORMS 

So far we have discussed only PN inspiral waveforms. 
To produce complete waveforms that include the late inspi- 
ral, merger and ringdown, we need to include results from 
numerical-relativity (NR) simulations. In this section we will 
show how the quadrupole-alignment procedure simplifies the 
production of hybrid PN-NR waveforms. 

A variety of methods have been introduced to construct 
hybrid waveforms for non-precessing configurations fPJl [331— 
l59l . and see Ref. [54] for a unified summary of the methods 
in use. In all methods the PN and NR waveforms are aligned 
at some time, or over a time or frequency window, and then 
blended together. Such waveforms have been used to produce 
phenomenological waveform models lfT2l[l3ll56ll57ll60l , and 
are now also being used to test GW search and parameter es- 
timation tools 11541 . A number of studies have also been per- 
formed on the length requirements of NR waveforms in order 
to produce sufficiently physically accurate hybrids ll59l |6TI - 
l63l and these also include estimates of the influence of errors 
and ambiguities in the hybridization process on the physical 
fidelity of the final waveform. 




FIG. 6: Evolution of J (red) and L (blue) plotted on the unit sphere, where J is initially aligned with the direction (0, 0,1). The left panel shows 
the evolution of these two directions for a case of simple precession. The precession cone described by J is very small in comparison to the 
one described by L. The right panel shows the same characteristic directions for a case of transitional precession. In this case J clearly moves 
along the unit sphere away from its initial direction and separates from L during the transition. 
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FIG. 7: The panel shows the magnitudes of the (2,2) modes for 
the transitional precession case before (red) and after (blue) the 
quadrupole alignment was applied. The change of the direction of 
/ at / = 1.587 • 10 6 M is indicated by the vertical line. A strong mod- 
ulation is introduced into the original waveform at that time, which 
is completely removed after quadrupole alignment. 



The construction of hybrids for precessing-binary configu- 
rations is more complex: not only do the time and phase of 
the PN and NR waveforms have to be aligned, but to some ex- 
tent the orientations of the spins and orbital plane must agree 
as well. For the precessing-binary hybrids that were used in 
Ref. lfT2ll . the hybrid waveforms were constructed by match- 
ing the NR waveforms with PN waveforms computed from 
the same PN evolution that was employed to construct the 
initial data for the NR simulations. This technique ignores 
mismatches in the binary orientation and physical parame- 
ters due to the emission of junk radiation [64 , 65 1 and gauge 
changes l66l l67l in the early stages of an NR simulation, al- 
though these effects are expected to be small; see Ref. ifTSl 
for a detailed discussion of this point in the context of non- 
precessing-binary hybrids. 
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FIG. 8: The two panels show the two Euler angles /3 and y deter- 
mined by the quadrupole-alignment procedure for the transitional 
case. The time when the z-component of J changes sign is indicated 
by the vertical line. 



These complications can be avoided through the use of QA 
waveforms. The PN and NR waveforms, both converted to 
the QA frame, can now be aligned exactly as in the non- 
precessing cases. In order to reverse the QA process, it is also 
necessary to align the QA angles (j3 , y, e), but this is straight- 
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forward, as we show below. 

In the next section we will outline how we produce a QA 
hybrid for the precessing-binary waveform that we used in 
Ref. ifTTl. This corresponds to the first configuration discussed 
in Tab. |lj q = 3, %\ = 0, %2 = 0.75, and S ■ L = 0. Having 
produced the QA hybrid, we will examine where our non- 
precessing-binary mapping hypothesis breaks down as we ap- 
proach merger. That the hypothesis must break down is clear, 
because the spin of the final merged black hole will be influ- 
enced by the black-hole spins in a way that the orbital phase 
evolution is not, and the mass and spin of the final black 
hole will not be the same as that for the corresponding non- 
precessing inspiral configuration. 



A. Construction of QA hybrids 

A QA hybrid can be produced by making use of the same 
procedure as for a non-precessing-binary hybrid. We will 
briefly summarize the method that we used. 

We start with a PN and an NR waveform, each for the same 
physical configuration. The last requirement is achieved to 
good approximation by using results from the PN evolution to 
produce the initial parameters for the NR evolution. The PN 
and NR waveforms are then put into the QA frame by the pro- 
cedure described in Sec. |IIC| We will produce a hybrid of ^4, 
and note that, since the QA frame is non-inertial, we cannot 
produce 'P^' 4 by taking two time derivatives of h9 A . We must 
first produce the 'P^m modes from the original precessing- 
binary GW-strain modes, him, and apply the QA algorithm to 
«P 4 ,2m. 

We then choose a matching frequency (0,„, and locate the 
times tpN and tug when each waveform passes through that 
frequency. For our q = 3 configuration, co m = 0.07. We 
align the PN and NR frequencies around that time such that 
both (j> PN (t PN ) = ^nr^nr) and (0 PN (t PN ) = (Onr^nr) = ®m- 
The hybrid waveform is then produced by blending together 

^4 A pN an d ^4 "nr w ith a linear transition function of width 
At — 200M around the matching frequency. The final wave- 
form is then 



jQA 



(4.1) 



where a± = (At /2 ± t ) / At when t 6 [—At, At] and zero or one 
otherwise, and the time has been shifted such that t = co- 
incides with the point at which co = CO,,,. This constitutes the 



QA hybrid. Fig. 10 shows the real part of ^4 around the time 
where the matching was performed, which is at t = 0. The 
figure shows the PN and NR waveforms, as well as the final 
hybrid, and we see that the matching between the PN and NR 
waveforms is smooth. 

To convert this into a physical precessing-binary 
hybrid, we also require hybrids of the QA angles 
{/3(f), /(?),£(?)}. These are produced as follows. 
The two angles {j5(t),y(t)} define a vector n{t) = 

{sin(-j3(0)cos(-y(0),sin(-j8(?))sin(-y( ? )),cos(-j3(0)}, 
which traces out a path on the unit sphere. The QA angles 
for the PN waveform define «/>^(f), while those for the NR 
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FIG. 9: Hybridization of the QA angles |3 and 7. Upper two pan- 
els: the black (dotted) lines indicate the inspiral PN values, the red 
(dashed) lines indicate the later NR values, and the green (solid) lines 
indicate the hybrids. The lower figure shows the evolution of the 
aligned QA directions, where here the black line indicates long PN 
inspiral of duration 2.9 x 10 5 M, and the red line indicates the NR 
results up to merger. 



waveform define nfjp(t). We perform a fixed rotation R/w to 
npN{t) (and another to n^ff(f)), such that both vectors 
are equal at the matching frequency, np^{tpN) — WAtR^W?)' 
Since the angle 7 is ill-defined when n = {0,0,1}, we do 
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300 



FIG. 10: The PN (red), NR (green) and hybrid (dashed black) wave- 
forms near the matching time (t = 0). The PN and NR waveforms 
are blended together in the window At = [—100, 100], indicated by 
the shaded region. 



not choose that as our (arbitrary) matching direction, but 
rather the vector such that j3(f/w) = 0.1 rad. Specification of 
a third Euler angle allows us to require that the vectors not 
only meet at the matching time, but that the curves they trace 
out are parallel at that time. To do this we simply measure 
the angle between the two curves at the matching time, and 
then rotate Hjvs(0 around the axis defined by the matching 
direction, njvsfevx)- Fig- [9] shows the first two angles at 
the times close to the matching frequency, and the final 
PN and NR curves are shown in the lower panel of Fig. [9] 
The hybrid angles are constructed by smoothly blending 
between the PN and NR angles, in the same way as for the 
QA waveform. The precessing-binary hybrid can then be 
constructing by performing the reverse QA procedure with 
{ 7 ,/M}^{-£,-/3,-7}. 



B. Breakdown of the non-precessing-binary equivalence 




FIG. 11: Matches between QA and non-precessing hybrids, for our 
standard q = 3 configuration. The horizontal axis represents the fre- 
quency at which both waveforms are cut off in the match calculation, 
and indicates that the two hybrids agree well (match > 0.97) right up 
to the merger, indicated by the vertical line. 



to Mco = 0.1, and is close to the point where PN waveforms 
are typically terminated in inspiral searches. Below this fre- 
quency the white-noise match is consistent with the results in 



Sec. Ill A The peak of the waveform occurs at fM = 0.07, 
which is indicated by a vertical line. The fiducial acceptable 
match of 0.97 is indicated by a horizontal line. We see that 
the match is at or above 0.97 through the merger, and only 
degrades significantly during the ringdown. 

Once again we emphasize that these matches were com- 
puted using a white-noise power spectrum. Nonetheless, these 
provide evidence that the QA procedure is valid very close to 
the merger, and perhaps even up to ringdown. We will discuss 
the implications of this result for waveform modeling in the 
final section. 



We expect the simple mapping between QA- and non- 
precessing-binary waveforms to break down near merger. As 
we have seen, the effect of the spins on the inspiral rate comes 
predominantly from the spin components parallel to the or- 
bital angular momentum; this is why our mapping works. At 
merger, however, the spin of the final black hole is, to first 
approximation, Jf m = L + Si +S2, where the orbital and spin 
angular momentum vectors are those at the point of merger. 
(A far more sophisticated treatment of the final spin ingredi- 
ents is given in Ref. [68 1, and a number of estimates of the 
final spin as a function of the initial configuration exist in the 
literature ||69Tl7T1l .) All components of the spin now become 
important and the appropriate parameterization may no longer 
be the effective total spin Xeff ■ 

It is instructive to investigate where the mapping breaks 
down, and we can use the hybrid waveform constructed in 



the previous section to do this. Fig. 1 1 shows the match be 



tween the QA hybrid constructed above, and a non-spinning 
q = 3 hybrid (which would be the appropriate non-precessing 
configuration during the inspiral). The match is calculated 
for a range of termination frequencies of the two waveforms. 
For reference, the frequency fM = 0.016 corresponds roughly 



V. DISCUSSION: A ROUTE TO GENERIC-BINARY 
WAVEFORM MODELS 

We have extended previous work on the quadrupole- 
alignment (QA) procedure to show that it can be used not only 
to cast precessing-binary waveforms in a simple form, but 
to map these waveforms onto a sub-family of non-precessing 
waveforms. We verified that this sub-family is parametrized 
by only mass ratio and an effective total spin parameter, and 
that the non-precessing waveform that best matches each QA 
waveform (with white-noise matches of at least 0.995), corre- 
sponds to our predicted Xeff value to within Ax < 0.04. The 
mapping was tested on a range of inspiral PN waveforms with 
mass ratios q = 3 and q = 10, and even on an example of 
transitional precession; in all cases the approximations holds 
well within the level of accuracy of the PN phasing. As a fi- 
nal test, we used the reverse QA procedure to "wrap up" a 
non-precessing-binary waveform, and found that it matched 
the corresponding true precessing waveform with a match 
of > 0.97 for all binary orientations. We also showed that 
this procedure can simplify the construction of hybrid PN-NR 
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waveforms, and that the approximate mapping seems to hold 
all the way through to merger. 

Our results suggest that generic precessing-binary wave- 
forms can be generated with good accuracy by applying the 
reverse of the quadrupole-alignment transformation to a small 
class of non-precessing-binary waveforms. These waveforms 
appear to faithfully represent the "true" precessing-binary 
waveforms up to the point of merger, and perhaps even up 
to the ringdown. The problem of constructing a generic wave- 
form model can then be factorized into the smaller problem of 
modeling the two rotation angles {j3 (f ), y(r )} as a function of 
the black-hole spins and the mass ratio. 

More concretely, we propose the following strategy: once 
the evolution of the Euler angles j3 (f ) and 7(f) has been deter- 
mined for a large sample of the configuration space, these can 
be modeled as functions that depend on some set of parame- 
ters X 

j8 = j8(X(0), (5.1) 

r = r(Mt))- (5.2) 

We emphasize that the X should be physical parameters, or a 
combination of physical parameters. The third angle e(t) is 
automatically determined given the two others. The rotation 
angles are unique up to an overall rotation of the frame of 
reference; we expect that they will assume the simplest form 
if in the limit of infinite binary separation = (0,0, 1). 

Since precessing inspiral-merger (IM) waveforms can be 
mapped onto non-precessing ones via quadrupole alignment, 
using the angles {y(f), /3(f), £(/)}, a phenomenological IM 
model with (anti-)aligned spins can be used and "twisted up" 
with the inverse angles {— e(f),— J3(f),— y(f)}. This will give 
us a phenomenological IM model, 

©)=RH -P -rX(i,*e ff ;0- (5-3) 

Needless to say, an inspiral model is not urgently needed: 
we can already produce generic waveforms by integrating the 
PN equations of motion, as we have done in this paper. A sim- 
ple closed-form approximation to these solutions could sig- 
nificantly improve the efficiency of gravitational-wave search 
and parameter-estimation pipelines, but there is no barrier in 
principle to producing theoretical estimates of any of these 
signals. The real need is for complete IMR waveform models. 

Given in addition a phenomenological model for the ring- 
down, hf m (Xj{;t), which is parameterized by some yet-to-be 

determined subset Xr of the full binary parameters X, we ex- 
pect that we can produce a combined IMR model, which can 
be indicated schematically as 

CW = R(-£-P-7)<(i,&;t)x4(i^). (5.4) 

For ease of use in GW searches, ideally this model would be 
cast in closed-form expressions in the frequency domain. 

We still have the problem of modeling a seven-dimensional 
parameter space, but we now have to model only two func- 
tions, and, as we can see from Fig.|3](and even Fig.[8]for tran- 
sitional precession), they are smooth, simple functions, that 



may be far easier to model than the complicated amplitude 
and phase modulations that are standard features of the physi- 
cal waveforms. It is also likely that many of the features of the 
full seven-dimensional parameter space can be captured by a 
model that considers only a subset of the parameters. It is also 
quite possible that we will need to employ a non-precessing 
model that treats both black-hole spins, and/or the effective 
spin that proves most useful will differ from that presented 
here. Our main purpose is only to outline a general research 
program to develop generic waveform models based on QA 
waveforms. 

In the short term, a number of issues deserve further 
study. One is that the QA method, in all forms proposed 
to date lfTTll45l [46l. has only been studied in detail prior to 
merger. It is likely that in order to apply the QA method opti- 
mally through the ringdown, it will be advantageous to make 
use of spheroidal (rather than spherical) harmonics, but that 
remains to be seen. Another will be the appropriate blending 
between the inspiral-merger regime and the ringdown, which 
is likely to be parameterized by the final mass and spin, and 
probably also the mass ratio and a second effective spin pa- 
rameter IT721 . It would also be instructive to explore the ef- 
fectiveness of such signals for both GW searches and param- 
eter estimation across a wide volume of the binary parameter 
space. 

We consider all of these to constitute a promising strategy 
for the construction of approximate generic waveform models, 
and will pursue them further in future work. 
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